Nonequilibrium phase transition by directed Potts particles 
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We introduce an interface model with g-fold symmetry to study the nonequilibrium phase tran- 
sition (NPT) from an active to an inactive state at the bottom layer. In the model, q different 
species of particles are deposited or are evaporated according to a dynamic rule, which includes the 
interaction between neighboring particles within the same layer. The NPT is classified according to 
the number of species q. For q = 1 and 2, the NPT is characterized by directed percolation, and 
the directed Ising class, respectively. For q > 3, the NPT occurs at finite critical probability p c , 
and appears to be independent of q; the q — oo case is related to the Edwards- Wilkinson interface 
dynamics. 

PACS numbers:05.70Fh,05.70Jk,05.70Ln 
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Recently, the problems of phase transitions in nonequi- 
librium systems have attracted considerable interest in 
the physical literature Q. For example, the nonequilib- 
rium phase transition (NPT) from an active state to an 
inactive state becomes one of the central issues in the 
field of nonequilibrium dynamics Q . It was shown that 
the number of equivalent inactive (or absorbing) states 
characterizes universality classes for the NPT || . For ex- 
ample, NPT, occurring in the monomer-dimer model for 
the catalytic oxidation of CO Q, the contact process ||, 
etc, has one absorbing state, and belongs to the directed 
percolation (DP) universality class ||. Most of dynamic 
problems exhibiting the NPT with absorbing state, be- 
long to this universality class. Meanwhile, there are a few 
exceptions: when there exist two absorbing states in the 
dynamic process, the critical behavior near the thresh- 
old of NPT is distinctive from that of the DP class, and 
belongs to the directed Ising (DI) class ||, equivalent 
to the class of parity-conserving branching-annihilating 
random walks The stochastic models in the DI class 
include, for example, the probabilistic cellular automata 
model Q| , the interacting monomer-dimer model H , the 
modified Domany-Kinzel model |To| . On the other hand, 
the scenario that the NPT can be classsified according to 
the number of absorbing states is not complete since the 
NPT with more than two absorbing states is not known 
yet. In other words, the NPT in the directed g-state 
Potts (DPotts) class for q > 3 is not discovered yet. One 
reason may be from that for q > 3, active sites (kinks) 
generated at the boundary between different states do- 
mains appear more often than in the DI case. Thus, the 
absorbing state could not be reached for finite control pa- 
rameter. Accordingly, it would be interesting to search 
for the NPT belonging to the DPotts class for q > 3. 

Recently, the NPT has been considered in association 
with the roughening transition (RT) from a smooth to 
a rough phase. The NPT behavior appears at a par- 
ticular reference height of the interface. For example, 
for the monomer deposition-evaporation (m-DE) model 
introduced by Alon et al , where evaporation of par- 
ticle is not allowed on terrace, the reference height is the 
spontaneously selected bottom layer. The site where the 



interface touches the reference height, called vacant site, 
corresponds to the active site of the NPT. In the active 
phase of the NPT, the interface fluctuates close to the 
reference height, being in a binded state, in which the 
interface is smooth. On the other hand, in the inactive 
phase, the interface is detached from the reference height, 
being in a unbinded state, and the interface is rough. 
Accordingly, RT accompanied by the binding-unbinding 
transition of the interface is related to the NPT at the 
reference height, which is characterized according to the 
number of symmetric states in the dynamic rule [ |l2fl . For 
the interface models with two-fold symmetry in their dy- 
namic rule, which arc the extensions of the m-DE model 
by assigning a couple of species to particles for one model 
p2| , called the two-species model, and by enlarging the 
size of particles for the other model jp| , called the dimer 
DE (d-DE) model, the dynamics at the reference height 
near the threshold of RT behaves similarly to the DI dy- 
namics. 

Both models contain the suppression effect against 
generating active sites, and the critical threshold of RT 
is considerably lower compared with the monolayer ver- 
sion. For the two-species model, the critical threshold is 
p c w 0.4480 for the interface version, but p c w 0.7485 
for the monolayer version. In this Letter, we introduce 
an interface model with arbitrary g-fold symmetry, called 
the g-species model, and examine the NPT at the refer- 
ence height for the cases q — 3, 4, 5 and oo. The result is 
compared with the previous one for q = 1 and 2. In ad- 
dition, we study RT for each q. Interestingly, it is found 
that the NPT at the reference height occurs at finite 
deposition-probability for all cases, and their character- 
istics for q > 3 arc different from the cases q = 1 and 2, 
but appear to be independent of q as long as q > 3. In 
particular, the interface model for the case q = oo is re- 
duced to the restricted solid-on-solid (RSOS) model with 
deposition and evaporation processes. In this case, since 
RT occurs when the probabilities of deposition and evap- 
oration are equal, the interface dynamics belongs to the 
Edwards- Wilkinson (EW) interface dynamics |14|] . Ac- 
cordingly, it is found that the NPT of the g-species model 
for q > 3 in one dimension is related to the EW class. 



In the g-species model, g different species particles are 
deposited or evaporated on one dimensional substrate 
with periodic boundary condition. A site is first selected 
at random at which either deposition or evaporation of 
a particle is attempted with probability p and 1 — p, re- 
spectively. In the attempt of deposition, one species is 
selected among the q species with equal weight 1/q. De- 
position or evaporation are realized under two conditions 
described below. First, the RSOS condition is imposed 
such that the height difference between nearest neighbor- 
ing columns docs not exceed one. Thus the attempt is 
realized as long as the RSOS condition is satisfied even 
after deposition or evaporation event. Secondly, the in- 
teraction between neighboring particles within the same 
layer is considered: When a particle of a certain species 
(e.g. A-species) is deposited on a hollow between par- 
ticles, the deposition is not allowed when both of the 
neighboring particle on each side are of a common species 
(B-species), different from the dropping one (A-species). 
Meanwhile, a particle of a certain species (e.g. A-species) 
is not allowed for evaporation when it is sandwiched be- 
tween particles of the same species (A-species) within 
the same layer. However, a particle can deposit or evap- 
orate when two neighboring particles on each sides are 
of different species from one another, or one (or both) 
of the neighboring sites is (are) vacant. As q — > oo, the 
probability of forming three particles in a row with the 
same species is zero. Thus, the secondary restriction is 
meaningless. Hence, the model is reduced to a random 
deposition-evaporation model under the RSOS condition. 
On the other hand, the initial substrate is flat where its 
height is referred as h = 0. The dynamics proceeds only 
for h > 0, so that evaporation of particles at h — is not 
allowed. 

Monte Carlo simulations are performed by varying the 
deposition probability p and system size L for the cases 
q = 3, 4, 5 and oo. We first discuss the dynamics oc- 
curring at the reference height, h = 0. We consider the 
density of the vacant site p(p, t) at the reference height 
with varying time t at a certain deposition probability p. 
When p is small, particles form small-sized islands which 
disappear after their short lifetime and the growth veloc- 
ity of interface is zero. The density p(p, t) saturates at 
finite value as t — > oo, as shown in Fig.l. As p increases, 
deposition increases and islands grow, until, above a crit- 
ical value p c , islands merge and fill new layers completely, 
giving the interface a finite growth velocity. Accordingly, 
RT occurs at p c , and the NPT at the reference height 
occurs as well. The critical probability p c , estimated for 
each q, is listed in Table 1. In particular, for q = oo, the 
critical probability p c is determined as p c = 0.5, when the 
probabilities of deposition and evaporation are equal. For 
p < p Cl the saturated value p(p) behaves as 

P(JP)~(PC-PY, (1) 

with the order parameter exponent f3. We estimated 
the exponent j3 from the data obtained from system size 



L = 500, which is listed in Table 1. The numerical values 
for q > 3 do not seem to be close to each other, which 
makes it hard to conclude that the cases of q > 3 belong 
to the same universality class. However, the exponent 
is hard to be measured precisely, because it is extremely 
sensitive to the estimated value p c . Accordingly, the nu- 
merical value of /? even for q = 2 is rather broadly ranged 
as can be noticed in Ref. (lOj. On the other hand, for 
p > Pc, the density decreases to zero exponentially in the 
long time limit and at the critical threshold p Cl it decays 
algebraically as 

p(p c ,i)~ (2) 

The power, (3/v\\, is measured in Fig.l and is tabulated 
in Table 1. The values are reasonably close to one an- 
other for q > 3, but considerably larger than the values 
for q = 1 and 2. This result suggests that the number of 
species is unimportant for q > q c = 3. On the other hand, 
we need to discuss the relation between the vacant sites 
and the kinks which are the domain boundaries between 
different species, because the kinks are indeed active sites 
in the directed Potts dynamics. We measure the kink 
density pK{pc,t) by counting the sites with height equal 
to one and whose neighbors are different species, which 
behaves similarly to p(p c , t), as shown in Fig. 2. Thus we 
confirm that the vacant site density is equivalent to the 
active site density in the g-species model. Finally, we 
simulate the monolayer version of the g-species model, 
and find that the system does not exhibit NPT and is 
always in the active state for finite p. A similar behavior 
was found in the three-species monomer-monomer reac- 
tion model introduced by Bassler and Browne fl5| , in 
which the system is always in a reactive phase when the 
absorption rates for each species are identical. 

The density p{p, t) can be thought as the return prob- 
ability Po(p,t) of interface height h(x,t) to its initial 
height h(x, t) = after passing time t, averaged over 
the substrate position x. The subscript means that 
the time is measured from t — 0. In general, Po(p,t) 
is different from the first return probability Fq (p, t) of 
the interface height to its initial height, which decays as 
FpOp, t) ~ t~ 6 where 8 is called the persistent exponent 
fll6fl . For the EW interface under the given boundary 
condition at h = 0, we obtain that Fo(p c ,t) ~ <~ a75 , 
similar to Po(p Cl t). Hence, the power f3/v\\ for q = oo is 
related to the persistent exponent for the EW interface. 

Next, let us consider the size dependence of p(p c , t). To 
do so, we study the averaged density p^(p c ,t) over the 
active runs which contain at least one vacant site at the 
reference height. Since p(p c , t) decays as Eq.(2) at p c , the 
vacant sites disappear completely as t —> oo. There exists 
a characteristic time r such that every run is active up to 
the time r, while for t > r, it is active with probability 
less than one and proportional to p(p Cl t), leading to the 
saturation of p^ s '(p c ,t), as shown in Fig. 3. This fact im- 
plies that for t > t, when vacant sites exist, their number 
should be of order of unity. Once they are occupied, the 
sample falls into the inactive state, because other sites 
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are covered by particles on upper layer. Hence, the den- 
sity follows p(p c ,t) w 1/L. The characteristic time r 
depends on system size L as t ~ L 2 with the dynamic 
exponent z = vn/v±. According to the scaling theory, 
the saturated value, 



p( s \p c ) is written as 



P (S) (PC 



-P/v± 



(3) 



The power f}/v±_ is found to be » 1, as predicted. There- 
fore, the dynamic exponent z is obtained to be between 
1.32 and 1.41, which deviates from the value zew = 2 
for the interface dynamics. The origin of this discrep- 
ancy may be that the dynamics at the reference level is 
affected by the boundary at h = 0. 

Let us consider the dynamics of interface above the 
reference height. We examine the interface fluctuation 
width, W 2 (L,t) = iJ2 x h 2 (x,t) - (^E^OM)) 2 at 
p c . Contrary to the cases q — 1 and 2, W 2 exhibits 
the power-law behavior as shown in Fig. 4, 



W 2 (L,t) 



[t 2 ^ for 
\L 2 *, fori 



(4) 



where £ and x are the growth and the roughness expo- 
nents, respectively. The exponents C and \ f° r 9 > 3 are 
obtained numerically as C ~ 0.22, 0.23, 0.23 and 0.24, and 
X « 0.46, 0.48, 0.46 and 0.48 for q = 3, 4, 5 and oo, respec- 
tively. These values are close to the EW values, £ = 1/4 
and x = 1/2- We also examine the height-height corre- 



lation function, C 2 (r) =< (h(r) - h(0)) 2 >r 



r 2x in the 



long time limit. The exponent x' is consistent with X- 
For p > p c , the scaling of the interface width belongs to 
the Kardar-Parisi-Zhang universality class fT^ ]. 

We also generalize the d-DE model into the trimer case, 
called the t-DE model, where a trimer is deposited with 
probability p or evaporated with probability 1 — p ac- 
cording to a rule similar to the d-DE model. Note that 
the t-DE model is different from the model introduced 



by Stinchcombe et.al. 18 in the two aspects: the for- 
mer (latter) model is an interface (monolayer) model, 
and evaporation on terrace is prohibited (allowed). As 
shown in Fig. 5, there exists a critical deposition prob- 
ability p c such that for p < p c , the density of the va- 
cant sites at the bottom layer is saturated, whereas for 
P > Pc, it decays algebraically. Note that the behavior 
of the exponential-type decay for p > p c does not appear 
in the t-DE model. The exponent, ft/vn is obtained to 
be k 0.38, which is almost half of the value for the 3- 
species model, which is also the case for q = 2. On the 
other hand, the dynamic exponent z w 2.47 is obtained 
by measuring the characteristic time for finite-size cutoff. 
This value is inconsistent with the one for the g-species 
model. Accordingly, further study for the dynamic ex- 
ponent z is required. Detailed numerical results for the 
t-DE model and its generalization into arbitrary g-mer 



case will be published elsewhere |19 



by the binding-unbinding transition with respect to the 
reference height. The NPT occurring at the reference 
height has been examined numerically for q = 3, 4, 5 and 
oo. We found that the NPT occurs at finite p c for all 
the cases of q, which is remarkable, because the NPT for 
q > 3 has never been found before. For q > q c = 3, 
the number of species is unimportant, and the NPT in- 
dependent of q is related to the EW interface dynamics 
via the return probability of interface to its initial height. 
We measured the critical exponents for NPT and RT. We 
also considered the t-DE model with three-state symme- 
try, and compared the result with the three-state model. 

The authors wish to thank H. Park and J.D. Noh for 
helpful discussions. This work was supported in part 
by the Korea Research Foundation (98-00 1-D00280 & 
98-015-D00090). 



In summary, we have introduced the g-species model, 
which is an interface model exhibiting RT accompanied 
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FIG. 1. Double logarithmic plot (DLP) of p(p,t) versus 
time t (a) for q = 3 at p = 0.4904, 0.4910, 0.4913, 0.4914(= p c ) 
and 0.4917 from the top. (b) for q = 4 at p = 0.5020, 
0.5024, 0.5026, 0.5027(= p c ), and 0.5031 from the top. (c) 
for q = 5 at p = 0.5051, 0.5059, 0.5063, 0.5065(= p c ), and 
0.5066 from the top. (d) for q = oo at p = 0.4993, 0.4997, 
0.4999, 0.5000(= p c ), and 0.5005 from the top. All data are 
obtained for L — 500, averaged over 1000 realizations. Insets: 
DLP of p(p) versus e for each q. The solid line is a guideline 
to the eye. 
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FIG. 4. DLP of W 2 versus time f at p c for different system 
sizes L = 100, 200, 300, 500, 800 and 1000 from the bottom for 
(a) q = 3, (b) q = 4, (c) q = 5, and (d) q = oo. Insets-left: 
DLP of C" 2 (r) versus distance r for L = 500 for each 5. The 
dashed line is a guideline to the eye. Insets-right: DLP of the 
saturated value of W versus L for each q. The solid line is 
a guideline to the eye. All the data are averaged over 1000 
realizations 
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FIG. 2. DLP of p(p c , t) (top) and the kink density pK{p c ,t) 
versus time t at p c . The numerical results are obtained for 
L = 1000, averaged over 1000 realizations. 




FIG. 3. DLP of p ( - s \p c ,t) versus time t at p c for different 
system sizes L = 100, 200, 300, 500, 800 and 1000 from the 
top for (a) q = 3, (b) q = 4, (c) q = 5, and (d) q = 00. 
The data are averaged over 1000 realizations. Insets: DLP of 
p' s '(pc) versus L for each q. The solid line is a guideline to 
the eye. 
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FIG. 5. (a) DLP of p(p, t) versus time t for the t-DE model 
at p = 0.32, 0.34, 0.36, 0.38(= p c ), and 0.42 from the top. The 
data are obtained from system size L — 500, averaged over 
500 realizations, (b) DLP of p(p c , i) versus time t for various 
system sizes L = 99, 198, 300, 498 and 798 from the left. The 
data are averaged over 500 realizations. Inset: DLP of the 
characteristic time r versus system size L. The solid line is a 
guideline to the eye. 



TABLE I. The estimated values of p c , /?, (3/v\\, f3/v±, and 
2 for various q. The data for the cases of q — 1 and q = 2 arc 
quoted from Refs.12 and 13, respectively. 



q 


Pc 


P 






z 


1 


0.189 


0.28 


0.16 


0.25 


1.56 


2 


0.4480 


0.88 


0.59 


0.97 


1.64 


3 


0.4914 


1.01(3) 


0.75(1) 


0.99(1) 


1.32 


4 


0.5027 


0.94(2) 


0.73(1) 


0.99(1) 


1.36 


5 


0.5065 


0.91(3) 


0.70(1) 


0.99(1) 


1.41 


00 


0.5000 


0.95(3) 


0.75(1) 


0.99(1) 


1.32 
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